Generalized equation of state for cold superfluid neutron stars 
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Abstract 

Mature neutron stars are expected to contain various kinds of superfluids in their interiors. 
Modeling such stars requires the knowledge of the mutual entrainment couplings between the 
different condensates. We present a unified equation of state describing the different regions of a 
neutron star with superfluid neutrons and superconducting protons in its core. 



PACS numbers: 97.60.Jd,26.60.-c,26.60.Kp 



The existence of superfluids inside neutron stars is rather well-established. It was first 

n 

suggested a long time ago [1[ before the actual discovery of pulsars and was later supported 
by the observation of the long relaxation time following the first observed glitch in the Vela 
pulsar j^J. Since then, glitches have been observed in many other pulsars. Superfluidity 
is also expected to play a key role in the origin of large glitches and several scenarios 
have been proposed, e.g. a eatastrophie unpinning of superfluid vortices g a transition 
between turbulent and laminar superfluidity || or a superfluid two-stream instability |5j. 
Microscopic studies [6j indicate that the interior of mature neutron stars contain three kinds 
of condensates: a neutron superfluid permeating the inner layers of the crust and a mixture 
of a neutron superfluid with a proton superconductor in the core. The most massive neutron 
stars may contain additional superfluid species in their central cores such as hyperons and 
quarks. The critical temperatures below which these various transitions occur still remain 
very uncertain. 

At finite temperatures, a superfluid is described by two velocity fields: one associated 
with the superflow and a normal fluid which carries all the entropy. This two-fluid model 
was first introduced by Tisza 7|] for explaining the unusual properties of superfluid helium 
4 He. With the discovery of superfluid 3 He, it was realized that in superfluid mixtures such 
as 3 He- 4 He, the different condensates are coupled by (non-dissipative) mutual entrainment 
effects sj]: even though they can flow with their own velocity field, the momentum of each 
superfluid is not aligned with its corresponding velocity. Entrainment effects may have an 
important impact on the evolution of superfluid neutron stars (see e.g. Ref. |9|). B. Carter 



has developed an action principle for obtaining the relativistic hydrodynamic 



equations 



of superfluids and superconductors, as found in the interior of neutron stars [10|. The 
superfluid mixture is described by a Lagrangian density A which depends on the 4-current 
vectors n£ of the different constituents labeled by X (/i,z/ are used to denote spacetime 
indices). Considering variations of the fluid particle trajectories, the action principle leads 
to the hydrodynamic equations 

« x Hunvx = /J» (i) 

where /* is 4-force density covector acting on the fluid X and the vorticity 2-form vu* 
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V^7Ty — V^7r^ is defined as the exterior derivative of the 4-momentum covector 

dA 



* = = !j t M^" ! ; + ^ Y K) ■ (2) 



rX 



The "anomalous" coefficients A XY above arise from entrainment effects. This variational 
ism has been adapted to non-relativistic superfluids using a 4D fully covariant frame- 



forma 

work JjJ in order to facilitate the matching between macroscopic and microscopic models, 
and also because Newtonian models have been widely used for qualitative studies of super- 
fluid neutron stars. The relation between relativistic and non-relativistic models has been 



discussed in Ref. 



12|. 



The electrically charged particles inside neutron stars are locked together by the interior 



magnetic field and co-rotate on very long time scales of the order of the age of the star [13 ]. 
In contrast the neutron superfluid being electrically uncharged can rotate at a different 
rate. This naturally leads to considering the interior of old neutron stars as a two-fluid 



mixture. However, as mentioned earlier 



may require more elaborate models [l4]. The simplest model of cold superfluid neutron 
stars thus consists of a mixture of a neutron superfluid with a plasma of charged parti- 
cles (superconducting protons and leptons in the core, nuclei and electrons in the crust). 
The core is assumed to be entirely superfluid and elasticity effects of the solid crust are 
neglected. The neutron and "proton" fluids are described by their 4-current vectors n£ and 
n% respectively. Introducing the particle densities n 2 c 2 = —g^n^n 1 ^ and using the notation 
x 2 c 2 = —g^n^rip, the Lagrangian density of the two fluids can be expressed to lowest order 



;he description of young and massive neutron stars 



in the relative currents as 



12| 



M n n> n p) = A (n n ,n p ) + X 1 (n n ,n p )(x 2 - n n n p ) . (3) 

The first term Ao = —S is related to the internal energy density £ in the absence of relative 
currents, while the second term accounts for entrainment effects. 

We have calculated the functions Xo(n n ,n p ) and Xi(n n ,n p ) for all regions of a neutron 



star using the self-consistent mean- field method with Skyrme effective interactions 15| . This 
method has been very successful in describing the structure and the dynamics of medium- 
mass and heavy nuclei and has been also widely applied to the description of neutron stars 
and supernova cores. This method allows for a unified treatment of both homogeneous 
and inhomogeneous matter with a reduced computational cost. In this work, we have used 



the Skyrme interaction BSkl7 underlying our Hartree-Fock-Bogoliubov (HFB) nuclear mass 



model HFB- 17 



16j |. This model fits essentially all the available experimental data with 



rms deviations of 0.58 MeV thus ensuring that nuclei in the crust of a neutron star will be 
properly described. Moreover this model was constrained to reproduce various properties of 
infinite homogeneous nuclear matter as obtained from many-body calculations using realistic 
nucleon-nucleon potentials. In particular, this model fits a realistic equation of state of 
neutron matter and can therefore be also reliably applied to describe the liquid core of 
neutron stars. We have determined the equilibrium structure and equation of state of the 
outer crust of neutron stars following the standard approach of Ref. [l^]. In this region, 
the only microscopic inputs are nuclear masses. We have used the HFB-17 nuclear mass 
table or experimental data when available. Results can be found in Ref. [IS]. For the inner 
crust, we have applied the fourth-order Extended Thomas- Fermi method with proton shell 
effects added via the Strutinsky-Integral theorem [l9| . We have found that above the density 
p ~ 1.4 x 10 14 g cm -3 , the crust dissoves into a uniform plasma of neutrons, protons and 
electrons. Muons appear at densities above p ~ 2.1 x 10 14 g cm~ 3 . Leptons are treated 
as relativistic Fermi gases. The functions X (n n ,n p ) and Xi(n n ,n p ) in the ground-state of 
cold dense matter are shown in Fig. „ (analytic expressions for the core can be found in 
Ref. [l2j). The mass-radius relation of non-rotating neutron stars (obtained after solving 
the Tolman-Oppenheimer-Volkoff equations) and the composition of their core are shown in 
Fig. q Causality is satisfied inside any stable neutron star. 

Even though entrainment effects in neutron star cores have been studied for a long 



time 
crust 



20 



it has been only recently realized that similar effects should also occur in the inner 



2l| . Taking into account the neutron superfluid in the crust is of prime importance for 



modeling pulsar glitches 



22] or quasi-periodic oscillations in Soft-Gamma Repeaters 23]. 



Calculations of Ai in the crust region will be reported elsewhere. 
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FIG. 1: Generalized equation of state of superffuid neutron stars, based on the HFB-17 mass 
model: binding energy per nucleon E/A = — Ao/^b — M where M is the nucleon mass (left panel) 
and entrainment coefficient Ai (right panel) as a function of the baryon density n^. 
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FIG. 2: Left panel: Mass-radius relation for non-rotating neutron stars using the unified equation 
of state shown in Fig. „ Right panel: fractions n x /nfe of protons (X = p), electrons (X = e) and 
muons (X = /i). 
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